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Abstract 

Recent angle-resolved photoemission (ARPES) experiments have indicated 
that the electronic dispersion in some of the cuprates possesses an extended 
saddle point near the Fermi level which gives rise to a density of states that 
diverges like a power law instead of the weaker logarithmic divergence usually 
considered. We investigate whether this strong singularity can give rise to high 
transition temperatures by computing the critical temperature T c and isotope 
effect coefficient a within a strong-coupling Eliashberg theory which accounts 
for the full energy variation of the density of states. Using band structures 
extracted from ARPES measurements, we demonstrate that, while the weak- 
coupling solutions suggest a strong influence of the strength of the Van Hove 
singularity on T c and a, strong-coupling solutions show less sensitivity to the 
singularity strength and do not support the hypothesis that band structure 
effects alone can account for either the large T c 's or the different T c 's within the 
copper oxide family. This conclusion is supported when our results are plotted 
as a function of the physically relevant self-consistent coupling constant, which 
show universal behavior at very strong coupling. 

PACS numbers: 74.20.Mn, 74.25. Jb, 74.20.Fg 
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On general grounds, we know that low-dimensional features in the electronic dispersion 
lead to singularities in the density of states, and we also expect that the superconducting crit- 
ical temperature T c will be enhanced when the Fermi level lies near one of these singularities. 
Early work along these lines focused on the A15 superconductors, where the low-dimensional 
features are one-dimensional chains, and calculations in both BCS0I and Eliashbergi~i for- 
malisms demonstrated a small but significant enhancement of T c . Current interest is in 
the family of high-temperature cuprate superconductors, where the low-dimensional struc- 
tures are the Cu0 2 planes. Energy dispersions from photoemission measurements and LDA 
band structure calculations led a number of workers to suggest that the density of states 
in these compounds was logarithmically divergent and that this divergence, in addition to 
a fairly typical pairing boson, could account for the large critical temperatures, the anoma- 
lous isotope effect, and the specific heat.H This picture of the origin of high-temperature 
superconductivity has become known as the "Van Hove scenario." While a logarithmic 
singularity is appropriate for La 2 _ x Sr x CuO^ (LSCO), recent photoemission measurements 
on YBa 2 Cu 3 6 . 9 (YBCO), YBa 2 Cu 4 8 , and Bi 2 Sr 2 CaCu 2 8 (BSCCO) suggest that the 
density of states in these materials has a much stronger power law divergence. 0~E1 Weak- 
coupling calculations suggest that the enhancement of T c due to these power law singularities 
would be much larger than that for the logarithmic singularity and may account for the larger 
T c 's observed in these compounds.0 

In this paper, we examine this scenario by using tight-binding fits to photoemission 
data as input to a strong-coupling Eliashberg T c calculation in order to determine whether 
power law Van Hove singularities significantly enhance T c 's over those given by logarithmic 
ones. We also quantify the T c enhancement as a function of singularity strength relative 
to a flat density of states. We find that the dramatic increase in T c due to the Van Hove 
singularity in weak coupling theory is moderated considerably by strong-coupling effects and 
that the enhancement of T c by a power law singularity over a logarithmic one is also reduced 
somewhat by strong-coupling effects. We suggest several possible extensions of this work, 
though, which should be investigated before drawing a negative conclusion concerning the 
Van Hove scenario for the cuprates. 

For computing the critical temperatures, we use the standard mean field formalismlll 
in which the electron self-energy is solved self-consistently from the single-exchange graph 
generalized to include the energy dependence of the density of statesBiH For simplicity, we 
take the pairing interaction to be an Einstein phonon with frequency Qq = 500 K (43 meV), 
and we assume that the electron-phonon matrix element g is independent of wave vector. 
The equations for the electron self-energy in Matsubara space H(iu n ) can then be written 
in the Nambu matrix notation^! as 



where G(E,iuj n ) is the electron Green's function which satisfies the Dyson equation 



with Go(E,iou n ) being the bare propagator. In these expressions, T is the temperature, g 2 
is the electron-phonon coupling constant, D is the phonon propagator, N(E) is the single- 
spin density of states per unit cell, r 3 is the third Pauli matrix, and uj n = (2n + 1)ttT are 
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Fermionic Matsubara frequencies. Throughout this paper, we set % = ks = 1- We solve 
these self-energy equations for all components of the Nambu self-energy (r , 7~i, and r 3 )ll! 
and at fixed band filling n determined by 



n 



]T e-^» 0+ J dE N(E) Tr [r 3 G(E, iu n )) . (3) 



We fix n in order to properly account for the loss of particle-hole symmetry in the interacting 
system, but we note that, due to the uncertainties in the parameterization of the band 
structures from ARPES data described below, the absolute values of n should not be taken 
too seriously. Note that this approach includes the effects of inelastic pair-breaking, which 
is important for obtaining reasonable T c estimates. We linearize these equations in the gap 
function and solve the resulting eigenvalue problem for T c . We also extract isotope effect 
coefficients a from the relation 



at — — 



1 d\nT c 



2 d\nQ 



(4) 



by changing the Einstein phonon frequency Qq a small amount and adjusting g 2 so that 
A = 2g 2 /Q W, W the band width, is constant. 

In order to complete the specification of the problem, we must discuss the band structures 
used to determine N(E) and the choice of coupling constants. We obtain the electronic 
dispersion E^ by fitting to a two-dimensional, tight-binding form E^ = Hf =0 ii?7i(k) with 
the basis functions r/i(k) listed in Tab. 1 and the parameters t{ in Tab. 2. Representative 
dispersions for LSCO and BSCCO are shown in Fig. 1. 

Our fitted band structures are obtained from different procedures, which we will now 
describe. For the LSCO case, a simple two-parameter fit was done to an LDA band structure 
calculation with a 2 eV band width. For the YBCO case, we perform three fits to the data. 
The first fit was done by assuming a 2 eV band width as in LSCO and adding a third tight- 
binding term to force a k 4 dependence along the k x direction around the (tc, 0) point. This 
procedure gives an extended saddle point and so approximates the photoemission data. 
A second fit was done by relaxing the third parameter to best fit all the data along (0, 0) 
to (tc, 0). This method splits the extended singularity at (tc, 0) into two nearby logarithmic 
singularities (bifurcated saddle points) as seen in LDA calculations, which are related to the 
displacement of the planar oxygens along the c axis of the crystaklll (The data are more 
consistent with an extended saddle point about (tc, 0) than bifurcated saddle points. The 
flatness of the experimental bands is almost certainly a correlation effect, which can only be 
approximately simulated by the tight-binding fit.) The third fit assumed a 1 eV band width 
and gives a better fit to the data along the (tc, 0) to (tc, tc) direction, which are consistent with 
a quadratic dispersion about (tc, 0); however, as this point is only 20 meV below the Fermi 
energy, little information is available before the band disperses through the Fermi energy. 
For BSCCO, the tight-binding fit was done to the proposed band structure presented by 
Dessau et al. based on their photoemission data and yields a band width of 1.4 eV. § The 
dispersion was fit so as to go like k A along both orthogonal directions about the (tc, 0) point. 
For comparison, we also include calculations for a flat density of states with a 2 eV band 
width. 
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In computing the density of states from these energy dispersions, we employ a two- 
dimensional "tetrahedron" code that produces data on an energy grid of spacing 0.1 meV 
near the singularity. We have checked that our results are insensitive to this mesh density 
and have reproduced our earlier calculations on a logarithmically divergent density of states 
where the sine; ularity is treated exactly.0 We find that, for the LSCO and bifurcated YBCO 
cases, the density of states has a logarithmic singularity. For the extended YBCO case, on 
the other hand, the divergence goes like the inverse fourth root of the energy difference from 
the singularity, and the divergence in BSCCO goes as the inverse square root. We remark 
that the densities of states are similar in both YBCO cases, even though the divergences 
are different. We see that, by studying these different materials, we are actually studying 
different Van Hove singularity strengths. 

The choice of coupling constant in our calculations is complicated due to the sharp 
structure in the density of states. In a material with a structureless density of states, the 
coupling strength is A = 2g 2 /WQ . However, the magnitude of the coupling extracted from 
transport measurements will be renormalized by singularities in the density of states, and 
so a more physical choice of coupling is to use the self-consistent ly calculated self-energy at 
the lowest Matsubara frequency ujq = nT c : Xz = — Im Ti{iujo) /uq. In what follows, we will 
show calculations with either A or Xz fixed and examine the differences. 

For completeness, we will discuss some caveats regarding our approach. First, our calcu- 
lations employ a mean-field theory where vertex corrections are ignored. Although justified 
by Migdal's TheoremS in conventional phonon-mediated superconductors with flat densities 
of states, this approximation may fail when the density of states is singular. Recent work 
has indicated that, for a logarithmic Van Hove singularity, vertex corrections may not be 
large and may actually reduce the computed T c ,0 but the case of a power law singularity has 
not been examined. Second, it has been pointed out that, because the electronic density of 
states enters into the screening of the electron-phonon matrix element, ignoring the detailed 
wave vector dependence of this matrix element and the phonon propagator-as we do-could 
lead to erroneous results. Third, phase-space restrictions on electron-electron scattering 
imposed by the Van Hove singularity enhances the electronic response functions, possibly 
leading to an electronic pairing mechanism.!!] Our calculations focus on a phonon-like pair- 
ing boson and so do not address this issue. Fourth, for YBCO, we consider only the plane 
band which contains the Van Hove singularity nearest the Fermi energy (for BSCCO, the 
splitting of the two CuO plane bands is ignored). We expect that this band will dominate 
the properties of the material, but full multiband calculations will be required to test this 
assertion. We therefore view the results of this paper as a first step in identifying the quan- 
titative effect of a singular density of states on T c when strong-coupling effects are included 
and leave these other extensions for future work. 

Having described our approach, we will now discuss our results, beginning with our weak- 
coupling calculations. For these computations, we drop the r and r 3 components of the self 
energy from the Eliashberg equations, and we fix A to 0.211 in order to obtain a 95 K T c for 
the YBCO band structure with the extended Van Hove singularity and 2 eV band width. 
Note that this weak-coupling calculation is not equivalent to the standard BCS calculation 
with a square well potential, but does serve to illustrate the effect of neglecting inelastic 
pair-breaking. 

The maximum T c occurs at a band filling near but not necessarily at the Van Hove singu- 
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larity, and we tabulate these maximum critical temperatures and isotope effect coefficients 
as a function of material type in Tab. 3. We observe that the maximum weak-coupling 
T c 's are enhanced by factors of 17 (LSCO) to 45 (BSCCO) over the flat density of states 
value and that this enhancement follows the strength of the singularity, increasing as the 
divergence goes from logarithmic for LSCO to fourth root for the extended YBCO to square 
root for BSCCO. In contrast, the isotope effect coefficient decreases systematically with 
singularity strength from the BCS result of 0.50 for the flat density of states to 0.20 for 
BSCCO. We also see that the extended and bifurcated densities of states for YBCO yield 
approximately equal T c 's, despite the differences in singular behavior. Reducing the band 
width by one half, though, reduces T c by 40%, despite working at fixed A. The isotope effect 
coefficient, however, is unaffected by this change. Finally, note that the critical temperature 
roughly doubles from the logarithmically divergent LSCO model to the fourth-root divergent 
extended YBCO model, similar to the T c 's of the actual materials. This observation leads 
to the speculation that the pairing interaction is the same in LSCO and YBCO with the 
difference in critical temperatures accounted for solely by the band structure. 

The strong-coupling calculations do not support this simple view, however. In Tab. 
3, we show strong-coupling results for the maximum T c obtained from the solution of the 
Eliashberg equations including all components of the self-energy and with A = 0.627, chosen 
once again in order to fix the T c for the W = 2 eV extended YBCO band structure to 95 K. 
We see that the critical temperatures still show a systematic increase with the strength of 
the singularity, but the strong-coupling T c 's show less sensitivity to the singularity strength 
than the weak-coupling T c 's: the maximum strong-coupling T c is enhanced by factors of 
only 2.6 (LSCO) to 3.4 (BSCCO) relative to the flat density of states T c . Similarly, the 
strong-coupling isotope effect coefficient still exhibits a systematic decrease with singularity 
strength, but the magnitude of this trend is severely reduced, leaving a roughly constant at 
0.4. We also note that the finite band width and strong-coupling effects combine to increase 
the flat density of states a to 0.51 instead of 0.50. As with the weak-coupling calculations, 
the extended and bifurcated densities of states yield similar T c 's, en's, and A^'s; reducing the 
band width to 1 eV reduces T c and Xz by approximately 20% but leaves a almost unchanged. 

Drawing these results together, we can say that, in both weak- and strong-coupling cal- 
culations, T c and a show systematic trends with the strength of the divergence in the density 
of states, but the magnitude of this effect is reduced in the strong-coupling calculations. In 
particular, we find that the isotope effect coefficient in strong coupling is almost indepen- 
dent of the strength of the singularity and is around 0.4, which is much larger than the 
observed en's in the optimally doped cuprates. Also, the enhancement of T c relative to the 
flat density of states in strong-coupling is two to three for our choice of A, as shown earlier 
for a logarithmically singular density of states.0 Furthermore, we find that these results are 
robust against the parameterization of the YBCO data (extended vs. bifurcated), but T c 
and Xz are found to be sensitive to band width (although a is not). 

Additionally, the idea of a universal pairing interaction does not seem to be supported by 
the strong-coupling calculations. By fixing the coupling in the extended parameterization 
of YBCO to give a strong-coupling T c of 95 K, we find that the maximum T c for LSCO 
becomes 78 K, much larger than the observed maximum T c of 40 K. This result casts doubt 
on the idea that the differences in T c are simply due to band structure effects. We should 
also remark that the experimental T c for BSCCO is lower than YBCO, despite the more 
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singular density of states inferred for the former. We show a plot of T c versus doping for 
this A in Fig. 2, which illustrates in greater detail both the reduced difference in T c and the 
enhancement of T c over that given by a flat density of states. 

There are several other features of Fig. 2 which should be pointed out. First, we 
note that the critical temperature rises around the van Hove singularity, and the width 
of this maximum broadens as the strength of the singularity increases. Away from the 
singularity, the T c falls below the flat density of states value. This behavior is a reflection of 
the redistribution of electronic states caused by the Van Hove singularity; basically, spectral 
weight is transferred from the band edge to the singularity, and T c follows this shift. Second, 
at low values of n, the plotted data show somewhat more complicated behavior due to the 
effect of the band edge. Finally, we see that the extended and bifurcated YBCO give similar 
T c vs. doping curves despite the differing singular behavior. The density of states peak in 
the extended case is higher than the bifurcated case but is also narrower in energy. This 
result supports the idea that what determines T c is an effective density of states which is a 
particular average of the actual density of states over the scale of the phonon energy f2 -§@ 

In looking at the self-consistent coupling strength Xz in Tab. 3, we see that Xz is larger 
than A for the singular densities of states by a factor of order two, implying that these 
calculations are in the strong-coupling regime. The systematic behavior of Xz with A is 
shown in Fig. 3. For a flat density of states and an infinite band width, Xz = X. For the 
finite band widths used here, this relation is valid only in the limit of small A. As seen in 
the inset to Fig. 3, at larger A, the flat density of states Xz is less than A by about a factor 
of two. On the other hand, when the density of states is singular, Xz is enhanced at small 
A, and the magnitude of the enhancement grows as A decreases and as the strength of the 
singularity increases. This behavior arises from the availability of more scattering states 
near the Van Hove singularity. At large A, however, the structure in the interacting density 
of states is smeared out by the strong scattering, and Xz for the singular density of states 
begins to behave more like that for a flat density of states; i.e., it is dominated by the finite 
band width. The curves come together at large A as shown in the inset to Fig. 3. 

Since both T c and Xz are enhanced by the Van Hove singularity, plotting the one quantity 
as a function of the other is instructive. To that end, we show in Fig. 4 the maximum transi- 
tion temperature T c maa; normalized by the Einstein phonon frequency f2 plotted against the 
self-consistent Xz for all the materials studied and with the results for a flat density of states 
included for comparison. Our results for the flat density of states reproduce the standard 
curve at small A^.0 Also at small Xz, T c is still strongly and systematically enhanced by an 
increasingly singular density of states compared to the results for the structureless density 
of states. This enhancement is reduced at stronger coupling, becoming at most around 50% 
for Xz ~ 1- At still stronger coupling, as seen in the inset to Fig. 4, the effect of the sin- 
gularity in the density of states is completely removed and the critical temperature follows 
the results for a structureless density of states. If T c ma:r /f2 is plotted against the bare A, 
the qualitative features of the curve are the same, but the T c enhancement is larger as one 
would expect from Tab. 3 and the results at very strong coupling do not scale as well as 
when the T c 's are plotted against Xz- 

We can understand these results in the following way. At small A, Xz is strongly enhanced 
by the presence of a singularity in the density of states, but it is approximately equal to 
A if there is no such structure. Comparison of critical temperatures at fixed Xz therefore 
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means that a much larger value of A is used in the flat density of states calculation and 
so the apparent effect of the singular density of states is reduced. At strong coupling, the 
enhanced scattering from the singular bare density of states acts to smear out the interacting 
density of states and so there is no T c enhancement. We note that the loss of the Van Hove 
singularity in the interacting density of states at strong coupling has been seen explicitly in 
recent calculations.il 

It is clear from our results that strong-coupling calculations are required to make physical 
predictions of the critical temperature. The remaining question is: should one compare 
these strong-coupling critical temperatures at fixed A or fixed A^? In our investigation of 
the possibility of a universal pairing interaction, it is the microscopic coupling g and phonon 
energy Qq which would be fixed, yielding-for a given band width-a fixed A. Our results in 
this case show little evidence for a universal pairing interaction due to the reduction of the 
effect of the Van Hove singularity on T c . Alternatively, one can extract an effective coupling 
constant from transport measurements which would include the self-consistent effects of the 
density of states and so could be used to fix Xz- From earlier work, EE we found that the 
canonical values of Xt r ~ 0.2 to 0.4 extracted from optical conductivity in YBCO correspond 
to Xz ~ 1, at least if spin fluctuations are the source of the resistivity. Similar estimates of Xz 
have appeared in the literature regarding phonon coupling strengths.^ Based on this result 
and Fig. 4, we conclude that, while the Van Hove singularity does enhance T c systematically 
with the strength of the divergence in the density of states, the enhancement factor in this 
parameter region is at most around 50% and so by itself cannot account for the large T c 's 
observed. In order to obtain a large T c , either a large pairing interaction energy Qq or a 
modification of the pairing interaction due to the singularity is required, but it not known 
whether either effect is consistent with the observed linear resistivity of the cuprates. 

In conclusion, we have performed Eliashberg computations of the critical temperature in 
superconductors with strongly singular densities of states chosen to model LSCO, YBCO, 
and BSCCO. We have demonstrated that, while the weak-coupling solutions do suggest a 
strong influence of the strength of the Van Hove singularity on T c and the existence of a 
universal pairing interaction in the cuprates, strong-coupling calculations show a reduced 
sensitivity to the singularity strength and do not support the notion that band structure 
effects alone can account for either the large T c 's or the different T c 's within the family 
of copper oxide superconductors. The reduced effect of the Van Hove singularities can be 
traced to the fact that the inelastic scattering included in strong-coupling calculations smears 
out the sharp structure in the density of states, producing a T c similar to that from a flat 
density of states. This effect is seen clearly when the results are plotted as a function of 
the physically relevant self-consistent coupling constant Xz and leads to universal behavior 
at very strong coupling. If the Van Hove singularities in the electronic band structures of 
the cuprates are responsible for high-temperature superconductivity, then the effect must 
come in through vertex corrections, screening effects in the electron-pairing boson matrix 
elements, direct electronic interactions, or multi-band effects. Based on the results of this 
work, investigation into these effects now acquires additional importance. 
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FIGURES 



FIG. 1. Electronic energy dispersion along high-symmetry directions in the Brillouin zone re- 
sulting from fits of (a) an LDA band structure calculation for LSCO and (b) ARPES measure- 
ments on BSCCO0 to the tight-binding form described in the text. The wave vectors in the 
two-dimensional Brillouin zone are T = (0,0), X = (tt,tt), and G (M) = (ir, 0). 

FIG. 2. Critical temperature T c in K as a function of the number of electrons per unit cell away 
from the optimum value n — n max computed for the tight-binding band structures corresponding 
to BSCCO (short-dashed line), bifurcated YBCO (long-dashed line), extended YBCO with 2 eV 
band width (solid line), and LSCO (dot-dashed line). The horizontal solid line is the value of T c for 
the flat density of states at half-filling. These results are produced by strong-coupling Eliashberg 
theory with an Einstein phonon of energy = 500 K and A = 0.627. See Tab. 3 for additional 
results at optimal doping n max . 

FIG. 3. Self-consistent coupling constant Xz as a function of the bare coupling constant A 
for the band structures corresponding to LSCO (solid boxes), YBCO with the extended (open 
diamonds) and bifurcated (open triangles) Van Hove singularities and 2 eV band widths, BSCCO 
(solid triangles), and a flat density of states of 2 eV band width (open boxes). The lines are to 
guide the eye. The A^'s result from the solution of the strong-coupling Eliashberg equations as 
described in the text. Inset: Data from the main figure displayed out to large A. 

FIG. 4. Maximum critical temperature T™ ax normalized by the Einstein phonon frequency fio 
as a function of the self-consistently determined Xz- Symbols are the same as in Fig. ^, as is the 
method of computation. Inset: Data from the main figure displayed out to large Xz in order to 
demonstrate the absence of density of states effects at very strong coupling. 
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TABLES 



TABLE I. Two-dimensional, tight-binding basis functions used in the fits of the energy disper- 
sion described in the text in units of the lattice spacing (a square unit cell is assumed). 



1 

1 \ (cos k x + cos k y ) 

2 cos k x cos ky 

3 \ (cos 2k x + cos 2k y ) 

4 \ (cos 2k x cos k y + cos k x cos 2k y ) 

5 cos 2k r cos 2k, , 



TABLE II. Parameters U in eV corresponding to the basis functions in Tab. 1 used to fit the 
LDA calculations and photoemission measurements of several cuprate superconductors. When ti 
is zero, this parameter was not included in the fit. For YBCO, in parenthesis are the band width 
(eV) and saddle point type (ext for extended, bf for bifurcated). 

LSCO YBCO(l,ext) YBCO(2,ext) YBCO(2,bf) BSCCO 



to 0.00 0.00 0.00 0.00 0.0607 

h -1.00 -0.50 -1.00 -1.00 -0.525 

t 2 0.20 0.15 0.38 0.38 0.100 

t 3 0.00 -0.05 -0.06 -0.09 0.0287 

U 0.00 0.00 0.00 0.00 -0.175 

t 5 0.00 0.00 0.00 0.00 0.0107 



TABLE III. Band filling corresponding to maximum critical temperature, n max , maximum 
critical temperature, T™ ax , isotope effect coefficient, a, and self-consistent coupling constant, Xz, 
for the densities of states discussed in the text. Results are obtained from weak-coupling Eliashberg 
theory with A = 0.211 and strong-coupling Eliashberg theory with A = 0.627 in order to give a 95 
K T c in the extended YBCO band structure with 2 eV band width. YBCO notation as in Tab. 2. 

Flat LSCO YBCO(l,ext) YBCO(2,ext) YBCO(2,bf) BSCCO 



n max 1.00 0.82 0.55 0.48 0.45 0.92 
Weak-coupling 

T rnax (jg 3 j 54 2 59.4 95.9 97.5 139 

a 0.50 0.34 0.28 0.28 0.27 0.20 
Strong-coupling 

T max 30 5 78 Q 77 9 95 Q 94 9 1Q5 

a 0.51 0.43 0.40 0.41 0.41 0.40 

X z 0.59 1.07 0.96 1.30 1.29 1.45 
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